On the origin of strong photon antibunching in weakly nonhnear photonic molecules 
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In a recent work [T. C. H. Liew and V. Savona, Phys. Rev. Lett. 104, 183601 (2010)] it was 
numerically shown that in a photonic 'molecule' consisting of two coupled cavities, near-resonant 
coherent excitation could give rise to strong photon antibunching with a surprisingly weak nonlinear- 
ity. Here, we show that a subtle quantum interference effect is responsible for the predicted efficient 
photon blockade effect. We analytically determine the optimal on-site nonlinearity and frequency 
detuning between the pump field and the cavity mode. We also highlight the limitations of the 
proposal and its potential applications in demonstration of strongly correlated photonic systems in 
arrays of weakly nonlinear cavities. 

PACS numbers: 42.50.Dv, 03.65. Ud, 42.25.Hz 



The photon blockade is a quantum optical effect pre- 
venting the resonant injection of more than one photon 
into a nonlinear cavity mode [l|, leading to antibunched 
(sub-Poissonian) single-photon statistics. Signatures of 
photon blockade have been observed by resonant laser 
excitation of an optical cavity containing either a single 
atom [2] or a single quantum dot [3| in the strong coupling 
regime. Arguably, the most convincing realization was 
based on a single atom coupled to a micro-toroidal cav- 
ity in the Purcell regime [J], suggesting that the strong 
coupling regime of cavity-QED need not be a require- 
ment. Concurrently, on the theory side there has been 
a number of proposals investigating strongly correlated 
photons in coupled cavity arrays [5|-l7| or one-dimensional 
optical waveguides [8|. The specific proposals based on 
the photon blockade effect include the fermionization of 
photons in one-dimensional cavity of arrays [9| , the crys- 
tallization of polaritons in coupled array of cavities [10| , 
and the quantum-optical Josephson interferometer in a 
coupled photonic mode system |11|. 



It is commonly believed that photon blockade necessar- 
ily requires a strong on-site nonlinearity U for a photonic 
mode, whose magnitude should well exceed the mode 
broadening 7. However, in a recent work |12l | Liew and 
Savona numerically showed that a strong antibunching 
can be obtained with a surprisingly weak nonlinearity 
(U ^ 7) in a system consisting of two coupled zero- 
dimensional (OD) photonic cavities (boxes), as shown in 
Fig. [TJ^a) [1^. Such a configuration can be obtained, 
e.g., by considering two modes in two photonic boxes 
coupled with a finite mode overlap due to leaky mir- 
rors: the corresponding tunnel strength will be desig- 
nated with J. In Ref. [12l | numerical evidence indicated 
that a nearly perfect antibunching can be achieved for 
an optimal value of the on-site repulsion energy U and 
for an optimal value of the detuning between the pump 
and mode frequency. However, a physical understanding 



of the mechanism leading to strong photon antibunching 
is needed to identify the limitations of the scheme in the 
context of proposed experiments on strongly correlated 
photons, as well as to determine the dependence of the 
optimal coupling and detuning on the relevant physical 
parameters J and 7. 

In this letter, we show analytically that the surprising 
antibunching effect is the result of a subtle destructive 
quantum interference effect which ensures that the prob- 
ability amplitude to have two photons in the driven cavity 
is zero. We show that the weak nonlinearity is required 
only for the auxiliary cavity that is not laser driven and 
whose output is not monitored, indicating that photon 
antibunching is obtained for a driven linear cavity that 
tunnel couples to a weakly nonlinear one. We determine 
the analytical expressions for the optimal coupling U and 
for the pump frequency detuning required to have a per- 
fect antibunching as a function of the mode coupling J 
and broadening 7. Our analytical results are in excellent 
agreement with fully numerical solutions of the master 
equation for the considered system. Before concluding, 
we discuss the experimental realization of such a scheme 
by using cavities embedding weakly coupled quantum 
dots. Moreover, we consider also the case of a ring of cou- 
pled photonic molecules showing that strong antibunch- 
ing persists in presence of intersite photonic correlations. 

We consider two photonic modes coupled with strength 
J; each mode has energy Ej and an on-site photon- 
photon interaction strength Ui [i = 1,2). The Hamil- 
tonian is written as 
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where hi is the annihilation operator of a photon in i- 
th mode, F and Wp arc the pumping strength and fre- 
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FIG. 1: (1) Sketch of the two coupled photonic modes. The 
couphng strength is J, and the antibunching is obtained with 
a smaU nonhnear energy U compared to mode broadening 7. 
(b) Equal-time second-order correlation functions g-. (r — 0) 
are plotted as functions of nonlinearity U — Ui = U2 normal- 
ized to 7. The nearly perfect antibunching is obtained at the 
pumped mode [g[f{T = 0) ~ 0] for U = O.O4287. Parame- 
ters: 71 = 72 = 7, J = 37, El — E2 = hujp + 0.2757, and 
Fi = O.OI7. 
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FIG. 2: (a) Optimal nonlinearity Uopt and detuning A_Eopt 
are plotted as functions of coupling strength J normalized to 
7 (71 = 72 = 7 and El = E2 — E). The perfect antibunching 
is obtained for J > 7/\/2- (b) Transition paths leading to the 
quantum interference responsible for the strong antibunching. 
One path is the direct excitation from |10) to |20), but it is 
forbidden by the interference with the other path drawn by 
dotted arrows. 



quency, respectively. Following Ref. [1^, we first cal- 
culate the second-order correlation function g\j {j) ~ 

{ix\aAT)aj{T)ai) / {a\ai){a,aj) in the steady state using 
the master equation in a basis of Fock states [13|. The 
results are shown as functions of nonlinearity U in 
Fig. [Ijb). As already demonstrated in Ref. [12|, we can 
get a strong antibunching of the pumped mode (gJi (0) ~ 
0) for an unexpectedly small nonlinearity U = O.O4287. 
In order to understand the origin of the strong anti- 
bunching, we use the Ansatz 



Coo|00)+c-'"p'(Cio|10) + Coi|01)) 

+ e-'2"p* (C20I2O) + Cnlll) + C02IO2)) + 



(2) 



to calculate the steady-state of the coupled cavity system. 
Here, \mn) represents the Fock state with m particles in 
mode 1 and n particles in mode 2. Under weak pumping 
conditions (Coo > Cw, Coi > C20, Cn, C02), we can cal- 
culate the coefficients Cmn iteratively. For one-particle 
states, the steady-state coefficients are determined by 



(A£;i - i7i/2)Cio + JCoi + FCoo 
(A£;2 - i72/2)Coi + JCio - 0, 



0, 



(3a) 
(3b) 



where A_Ej = Ej — hup and we consider a damping with 
rate 7^ in each mode. Since we assume weak pumping, 
the contribution from the higher states (C20, Cu, and 
C02) to the steady-state values of Cio,Coi is negligible. 
From Eq. (j3bp . the amplitude of mode 2 can be written 
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indicating that for strong photon tunneling (J 3> 
I Ai?2 1 , 72 ) , the probability of finding a photon in the aux- 
iliary cavity is much larger than the driven cavity. 



In the same manner, the coefficients of two-particle 
states are determined by 

2{AEi +Ui- i7i/2)C2o + V2JCn + V2FCw = 0, 

(5a) 

(A^i + AE2 - m/2 - i72/2)Cii + V2JC20 + V2JC02 

+ FCoi - 0, (5b) 

2{AE2 + U2- ij2/2)Co2 + V2JC11 = 0. (5c) 

When we simply consider Ei = E2 ~ E, and 71 = 72 = 
7, the conditions to satisfy C20 = are derived from 
Eqs. (gl) and ^ as 

j^{3AE + U2) - AAE^{AE + U2) = 2J^U2, (6a) 
12A£:^ + 8AEU2 - 7^ = 0. (6b) 

For fixed J and 7, from these equations, the optimal 
conditions (those that lead to C20 ~ 0) are given by 
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and, if J 3> 7, they are approximately written as 
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In Fig. [2ja), the optimal Ai^opt and Uopt [Eq. (O] are 
plotted as functions of J/7. The strong antibunching 
can be obtained even if U2 < 7, provided J > j/V2- 
Remarkably, the required nonlinearity decreases with in- 
creasing tunnel coupling J obeying Eq. (j8b|) . 
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FIG. 3: (a) The time evolution of the second-order correlation 
function, which oscillates with period 2n/J as the result of 
amplitude oscillation between |01) and |10). (b) Equal-time 
second-order correlation functions are plotted as functions of 
AEi = AE2 = AE normalized to 71 = 72 = 7. The spectral 
width of the antibunching resonance is « 0.87. Parameters: 
J = 37, Ui = U2= 0.04287, and F = O.OI7. AE = 0.2757 in 
panel (a). 



In Fig. mb), wc show a sketch of the quantum interfer- 
ence effect responsible for this counter-intuitive photon 
antibunching. The interference is between the following 

two paths: (a) the direct excitation from |10) — > |20) 
(solid arrow) and (b) tunnel-coupling-mediated transi- 
tion |10) A |01) A (111) O |02)) A |20) (dotted ar- 
rows) . In order to show in detail the origin of the quan- 
tum interference, we rewrite Eqs. (|6]) for C20 = as 
follows. First, we calculate Cn from Eqs. ^ and Q 
neglecting C20 as 



Cii = -2JFCio{AE + U2- ij/2){AE - i-f/2)-^ 

[2J2 - AAE{AE + U2) -f 7^ -h 127(2 A£: -f- f/2)]"^ 

(9) 

This amplitude is the result of excitation from |01) to 
|11) and of the coupling between |10) and |01) and also 
between |11) and |02). From this amplitude, C20 is de- 
termined by Eq. (|5ap as C20 oc JCn + FCiq, and we can 
derive Eqs. (O by the condition C20 = 0. 

As seen in Fig.[ljb), while no more than one photon is 
present in the first cavity mode at the optimal condition, 
there can be more than one photons in the whole system. 
While there is nearly perfect antibunching in the driven 

mode [ffii (t = 0) << 1], the cross-correlation between 

(I) 
the two modes exhibits bunching [§12 ij = 0) > 1]. The 

amplitude oscillation between |10) and |01) produces the 

time oscillation of 5ii (t) with period 2ii/ J as reported 

in Ref. [l^ and shown in Fig.[31[a). 

The equal-time correlation functions is plotted in 

Fig. [3]Jb) as a function of the pump detuning AEf-f. 

while the optimal value of the detuning is at AE = 

0.2757, a strong antibunching is obtained in a range of 

about 0.87 around the optimal value and the width of this 
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FIG. 4: Equal-time correlation functions are plotted as func- 
tions of coupling strength g between a cavity and a quantum 
dot. Solid line represents the results in the system sketched 
in the inset [Eq. (fTU]) ]. Parameters: 71 = 72 = 7cx ~ 7, 
J = 37, Si = fi^p + 0.2757, E2- E^ = 7, £ex - -Ba = 27 
and F = O.OI7. Dashed line represents the result in the sys- 
tem with one quantum dot and one cavity [Jaynes-Cummings 
model]. 71 = 7ex = 7, -Ecx — -Bi = 27, F = O.OI7, and fitjp 
is tuned to the lower one-particle eigenenergy of the Jaynes- 
Cummings ladder. 



window does not significantly depend on J/7. This may 
suggest that pump pulses of duration Atp longer than 
1/(0.87) could be enough to ensure strong antibunching. 
However, the timescale over which strong quantum cor- 
relations between the photons exist is on the order of 
1/J < ■\/2/7, as seen in Fig.[3lja). While weak nonlinear- 
ities do lead to strong quantum correlations, these corre- 
lations last for a timescale that scales with 1/J ex ^JUo■pt 
(see Eq. (|8b[) ). From a practical perspective, a principal 
difficulty with the observation of the photon antibunch- 
ing with weak nonlinearities is that it requires fast single- 
photon detectors IJ]. Conversely, for a given detection 



set-up, the required minimal value of the nonlincarity is 
ultimately determined by the time resolution of the avail- 
able single photon detector. 



As seen in Eq. ([6]), the nonlincarity Ui of the pumped 
cavity mode is not essential for the antibunching. This 
means that only the auxiliary (undriven) photonic mode 
must have a (weak) nonlincarity to achieve the quantum 
interference leading to perfect photon antibunching. As 
a practical realization, one could consider two coupled 
photonic crystal nanocavities, where the auxiliary cavity 
contains a single quantum dot that leads to the required 
weak nonlincarity (see the inset in Fig. 2]). The Hamil- 
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FIG. 5: (a) Sketch of a triangular lattice of coupled photonic 
'molecules'. The driven cavities (i = 1, 2, and 3) are coupled 
with strength J2 . (b) The equal-time second-order correlation 
functions in each mode (solid line) and between neighbors 
(dashed line) are plotted versus J2/7- Parameters: J = 87, 
Ei^hujp + O.45O7, U = 0.07697 and F = O.OI7. 



tonian is written as 
2 

-ffcav-JC = ^ Eiaja^ + J{ala2 + alai) 

+ ^ox|ex)(ex| + g (4lg> (ex| + H.. 

+ Fe-''^p*al + i^*e'"p*ai. (10) 

Here, |g) and |ex) represent the ground and excited states 
of the quantum dot, respectively, E^^ is the excitation 
energy, and g is the coupling energy vifith cavity mode 
2. Since the required nonlinearity is relatively weak, one 
can use a quantum dot which is off-resonant with respect 
to the cavity mode (|-Eox — ^2! > 72 =7) and/or does 
not satisfy strong coupling condition (t/ ~ 7). We take 
the quantum dot exciton broadening to be equal to the 
cavity decay rate for simplicity. We have solved numeri- 
cally the master equation associated to the Hamiltonian 
in Eq. (|T0| . Fig. |4] shows g\{{T = 0) of the pumped 
mode as a function of g/7. The coupling energy between 
the two cavities is J = 87, and then the required non- 
linear energy should be t/opt = 0.04287 from Fig. [2l In 
the present system, this nonlinear energy is practically 
achieved sX g = I.47, which is an intermediate strength 
between the weak- and strong-coupling regime of cavity 
mode and quantum dot excitation. The dashed line in 
Fig. |4] represents the results in the system consisting of 
one quantum dot and one cavity: in this ordinary Jaynes- 
Cummings system, only a small antibunching is obtained 
at 5' ~ 7, and the strong-coupling g ^ 7 is required for 
the observation of large photon antibunching [l|, |2[ . In 
contrast, in the new scheme using the quantum interfer- 
ence, a nearly perfect antibunching can be obtained even 
for g ~ 7. 

Finally, we note that the quantum interference can 
be generalized to a system of many coupled photonic 
molecules: in this case, the strong on-site antibunching 



can show an interesting interplay with quantum correla- 
tion between neighboring photonic modes. As a demon- 
stration, we consider a ring of three molecules whose 
driven dots are coupled with each other by a tunnel 
coupling of amplitude Ji [see Fig. [Slja)]. Also in this 
case a nearly perfect antibunching can be observed in 

(21 

each driven mode, as shown in the plots of g\^ (t = 0) 
as a function of J2/7 that are shown as a solid line 
in Fig. [5jb). In order to optimize the antibunching at 
a finite value of J2 ~ 7, values of C/ = 0.07697 ^'^^ 
AE = O.45O7 slightly different from the single-molecule 
optimal ones (t/opt = 0.04287 and AE'opt = 0.2757) had 
to be chosen. At the same time, a strong bunching ef- 
fect is observed in the equal-time cross-correlation func- 
tion between neighboring cavities, which shows a value of 

(2) 
glJAO) significantly larger than the coherent field value 

(2) 
of glJAO) = 1. This remarkable combination of strong 

on-site antibunching and strong inter-site bunching sug- 
gests that this system may be a viable alternative to the 
realization of a Tonks-Girardeau gas of fermionized pho- 
tons discussed in Ref. [9|. 

In summary, we have analytically determined that a 
destructive quantum interference mechanism is responsi- 
ble for strong antibunching in a system consisting of two 
coupled photonic modes with small nonlinearity {U < 7). 
The quantum interference effect occurs for an optimal on- 
site nonlinearity [/opt — tt^ jt , where J is the intermode 
tunnel coupling energy and 7 is the mode broadening. 
This robust quantum interference effect has the pecu- 
liar feature that the resulting quantum correlation be- 
tween the generated photons survive for timescales much 
shorter than the photon lifetime. Nonetheless, we have 
shown that this quantum interference scheme has the po- 
tential to generate strongly correlated photon states in 
arrays of weakly nonlinear cavities. 
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